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ABSTRACT 

We analyse the sample of pulsar proper motions, taking detailed account of the selec- 
tion effects of the original surveys. We treat censored data using survival statistics. 
From a comparison of our results with Monte Carlo simulations, we find that the mean 
birth speed of a pulsar is ~ 250 — 300 kms^^, rather than the 450 km s~^ found by 
Lyne & Lorimer (1994). The resultant distribution is consistent with a maxwellian 
with dispersion = 190kms~^. Despite the large birth velocities, we find that the 
pulsars with long characteristic ages show the asymmetric drift, indicating that they 
are dynamically old. These pulsars may result from the low velocity tail of the younger 
population, although modified by their origin in binaries and by evolution in the galac- 
tic potential. 
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1 INTRODUCTION 



The fact that pulsars have velocities much in excess of those 
of ordinary stars (a subset of whom are presumably the pul- 
sar progenitors) has been known almost since their discovery 
(Minkowski 1970; Trimble 1971; Lyne, Anderson & Salter 
1982). The origin of these velocities is not so clear. One pos- 
sibility is that they result from the disruption of a binary 
population (Gott, Gunn & Ostriker 1970; Iben & Tutukov 
1996), leaving the pulsar with a velocity characteristic of the 
orbital velocity of the progenitor in the binary. The prob- 
lem with this scenario is that it has trouble explaining the 
largest observed velocities (e.g. Phinney and Kulkarni 1994; 
see, however, Iben & Tutukov 1996). Another possibility is 
that the pulsar acquired its velocity from an asymmetric 
supernova collapse. This hypothesis has recently been bol- 
stered by the observations of binary pulsar PSR J0045-7319, 
in which the observed spin-orbit precession and orbital pe- 
riod decay are strong arguments for a natal kick (Lai, Bild- 
sten & Kaspi 1995; Lai 1996). 

Lyne & Lorimer (1994) have analysed the known sam- 
ple of pulsar velocities in the light of recent proper motion 
studies (Harrison, Lyne & Anderson 1993) as well as the 
new pulsar distance scale due to Taylor and Cordes (1993). 
They conclude that pulsars are born with a mean speed 
of ~ 450kms~^. Although Lyne and Lorimer restrict their 
sample to those younger than 4 x 10^ years to avoid the 
'vertical leakage' selection effect pointed out by Helfand & 
Tademaru (1977) and Cordes (1986), they did not treat the 
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selection effects that result from the flux limits of the pulsar 
surveys or the limiting accuracy of proper motion determi- 
nations. We shall attempt to do that here. Recently, Iben 
and Tutukov (1996) have also addressed this question, but 
in a less systematic fashion than what we propose to use in 
this paper. 

We restrict our analysis to only those pulsars with ve- 
locities determined by proper motion measurements. While 
scintillation data have been used to estimate velocities 
(Cordes 1986; Harrison and Lyne 1993) to within a factor 
of 2, we prefer to keep our sample as homogeneous as pos- 
sible, and so we exclude these data. We use the properties 
of well-known pulsar surveys to estimate the V^/Vmax correc- 
tion (Schmidt 1968) for each pulsar, making use of survival 
statistics (Feigelson & Nels on 1985) to treat those data with 
upper bounds (section 2.2). This allows us to estimate the 
two dimensional velocity distribution of the observed pulsar 
population and thus the kick velocity distribution taking 
into account the differential galactic rotation in section ^. 



2 THE PROPER MOTION DISTRIBUTION 
2.1 Selection Effects 

The Princeton pulsar database (Taylor, Manchester & Lyne 
1993) now contains ~ 800 pulsars, 101 of whom have mea- 
sured proper motions or upper limits. The number of surveys 
responsible for this profusion is also gradually increasing in 
size (in excess of 15). To do a proper treatment of the se- 
lection effects for the full proper motion sample would then 
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require modelling the selection effects of a significant num- 
ber of these surveys. Luckily, we note that most of the proper 
motion pulsars were detected in the earlier surveys. By re- 
stricting ourselves to those pulsars detected in the Molonglo 

2 (Manchester et al. 1978) and Green Bank/NRAO 1,2 and 

3 surveys (Damashek et al. 1978, Dewey et al. 1985, Stokes 
et al. 1985, Stokes et al. 1986), we are left with 86 out of 
101 pulsars. We note that 12 of the 15 pulsars left out have 
P< 0.1 s (and only one of our restricted sample satisfies this 
criterion), which is not surprising, since many of the later 
surveys focussed on finding faster spinning pulsars. In per- 
forming this cut, we lose one young pulsar and all but one of 
the pulsars with characteristic ages greater than 1 Gyr. This 
also means we are not affected by the possibly different evo- 
lutionary histories of millisecond pulsars (in particular the 
influence of binaries) . We shall omit the one old binary pul- 
sar (PSR 0655+64) which does fall into our sample as well. 
Figure |l| shows the distribution of the included and excluded 
pulsars as a function of period and velocity. 

The observed pulsar sample suffers from two obvious 
selection effects, due to flux limits and proper motion lim- 
its respectively. Figure |^ shows the distribution of inferred 
luminosities and transverse velocities. The lack of faint, fast 
pulsars (upper left corner) and bright, slow pulsars (lower 
right corner) is evident. 

To correct for this bias, we need to weight the pulsars 
according to the maximum volume in which they could have 
been detected, i.e. using a Vmax weighting. To do this we 
need to consider the detection efficiency of the various pulsar 
surveys. After Narayan (1987) (see also Dewey et al. 1984), 
the minimum flux detectable is 



(Tree ~h ^sky) 

To 



PW 



We (P-W)' 



(1) 



where Tree and Tsky are the system and sky noise tempera- 
tures, 5*0 is the flux normalization, and W is the measured 
pulse width. We is the intrinsic pulse width, which is broad- 
ened because of sampling, dispersion and scattering so that 



W'^ = We + r,in,p + Tom + ^s^catt- 



(2) 



Thus we calculate the minimum flux for a given survey in a 
particular direction. The parameters describing each survey 

were taken from Narayan (1987)^ and Stokes et al. (1986). 
We used the updated electron distribution model of Taylor 
and Cordes (1993) to calculate the dispersion and scatter 
broadening along a given line of sight. This accounts for the 
bias due to the flux limits. No such simple model exists for 
treating the proper motion limits. This is because the accu- 
racy of a given proper motion measurement depends on the 
vagaries of the distribution of background radio sources near 
the pulsar position on the sky (Harrison, Lyne & Anderson 
1993). As a crude model of this, we flt the distribution of 
proper motion errors in Harrison et al. using the distribution 
p(/i) = exp(— ^/10.5 masyr"^). In the Vmax calculations to 
follow, the proper motion cutoff is randomly selected from 
this distribution for each line of sight and . We also used 



t Narayan's formula for the sky background temperature in equa- 
tion (3.5) contains a typographical error. The factor + (b/3)'^j 
should be in the denominator. 



simple limits of 5 and 2 masyr"^ as a test of the impor- 
tance of this selection effect. This introduces a variation of 
~ lOkms"^ in the mean velocity of the young sample and 
about 30kms~^ in the old sample, indicating that most of 
the Vmax are flux-limited rather than proper motion lim- 
ited. Iben and Tutukov (1996) have also taken account of 
this selection effect. 

For each pulsar, we randomly place it in different direc- 
tions and at different distances with respect to the observer 
and calculate whether or not it would be detectable with 
it's known luminosity and transverse velocity (uncorrected 
for the local standard of rest) in any of the surveys we con- 
sider. Thus, using this Monte Carlo integration procedure, 
we determine the volume within which each pulsar could 
have been detected. These Vmax values determine the rela- 
tive weights of each of the pulsars in the corrected sample. 

A possible source of concern with this procedure is illus- 
trated by Figure ^. The analytic flux limits do not describe 
the complications of the true detection limits perfectly (in- 
terstellar scintillation makes the apparent flux vary). In Fig- 
ure ^ we find four pulsars detected by the Molonglo survey 
that lie below the analytically described detection thresh- 
old for that survey. This will reduce the weight accorded 
to these pulsars. To estimate the impact of this error on 
our results, we repeated the analysis with these pulsars ar- 
tificially 'brightened' to meet the flux limit expression. The 
mean proper motion we infer for the young pulsar sample 
increases by only 7 kms~^, so this is not a signiflcant source 
of error for our analysis here. 

The Vmax correction is not without biases of its own. 
In particular, weighting pulsars by their Vnax presupposes 
that the real population is distributed uniformly through- 
out the galactic volume. However, the pulsars are born from 
a disk population with a scale height of about 150-450 pc 
(Narayan and Ostriker 1990) . Thus, a population born with 
small velocities will not expand to fill as much of the spher- 
ical volume as a fast population. The above analysis then 
overcorrects for the slow pulsars (see Helfand and Tademaru 
1977; Cordes 1986; Lyne & Lorimer 1994). In order to adjust 
for this, we consider the maximum detectable volume to be 
limited in the vertical extent by the scale Zma,x=Vttp, where 

1/3 

tp is the pulsar timing age. If this is larger than Deq ~ Vmax, 
then there is no change in the weight assigned to that pulsar, 
but, if Deq > Zina.x, then we assume that we see the edge 
of the distribution of pulsars of this velocity, and reduce the 
weight given to that pulsar accordingly(see Figure ^ . Fur- 
thermore, for young pulsars, we set a lower limit on Zmax of 
450 pc, representative of the initial scale height. 

In Figure | we show the distribution of velocities with 
age and the relative weighting of each pulsar. 



2.2 Survival Statistics 

The Vmax correction takes care of the selection effects, but 
we still need to account properly for those data which only 
have upper limits ('censored' data in the statistical lexi- 
con). Of our 85 pulsars, 20 fall into this category. Using only 
those data with actual detections will bias our distribution 
to higher values as we will see below. 

Following Feigelson and Nelson (1985), we use survival 
statistics to treat the effects of our censored data. In partic- 
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ular, we use the Kaplan-Meier estimator (Kaplan and Meier 
1958) to calculate the cumulative probability distribution of 
the transverse velocities. However, we need to modify this 
method slightly to take account of our Vmax correction. 

Consider our data to be the set of n distinct values 
{xi}^^^, where the Xi are ordered in the manner x\ > X2 > 
■ ■ ■ > Xn, where Xi can be either a detected value or an upper 
bound. Now, for a given value t, let 

P^=P[t< x,+i \t<x^], 

i.e. the conditional probability that t is less than Xi+i given 
that it is less than Xi. 

Using this, we can calculate the probability that t is less 
than any given Xi {i > 1) by 

P[t< X,] = TT^jZlP [t < Xj+i \t <Xj]. 

To calculate this, we need to estimate the Pj. To start. Pi = 
1, since all values are at most as big as xi. For any other Xj, 
if it is a detection, then there are n — j + 1 values at most 
as large as Xj, and all except Xj are also at least as large as 
Xj+i. Thus, we estimate 

D _ " - i _ 1 1 



n ~ j + 1 n — j + 1 

If Xj is an upper bound, then Pj = 1 by following similar 
reasoning as above. Thus, we have that 

P[t<x,]=U'jZ\ f 1- 



where Sj = 1 if X j IS 3, detection, and S{j) = if Xj is 
an upper limit. This is the Kaplan-Meier (1958) estimator 
of the distribution function. In the case where we have ties 
in our data (i.e. more than one measurement at the same 
value), this becomes 



1 _ !k. 



where rij = number of measurements < Xj, and dj is the 
number of measurements at the value Xj. 

We also want to calculate the mean of velocities. Since 
p{x)dx = AP, the mean of a quantity x is given by 



<x>^ 



xdP - 



Pdx, 



where we have integrated by parts. Thus we estimate the 
mean by 



<Vt >=^P[x,] {x. 



and we take xq = 0. 

To include a V/Vmax weighting, we adjust the number 
of 'counts' at each value according to the weights w(j), and 
thus we have 

Uj ^-^ Wi 

i<j 

In calculating the error on this new estimator, we note that 
we have artificially increased the total number of observa- 
tions fed into the sum and thus have reduced the error by a 
factor \/Nw, where A'w is the number of pulsars including 



the weigh tings. W e remove this bias by multiplying the er- 
rors by y'Nw/N where is the true number of points in 
the sample. 



2.3 The Corrected Distribution 

Figure |^ shows the cumulative probability distribution of 
observed transverse velocities taking into account different 
levels of adjustment. With no adjustments, the mean trans- 
verse velocity of the entire sample is 355 km s~^, in agree- 
ment with the analysis of Lyne & Lorimer (1994). The mean 
of the entire sample including all corrections is 195 km s~^. 
However, if we restrict the sample to pulsars with charac- 
teristic ages < lO'^ years, we obtain a mean transverse ve- 
locity of 237 ± 49 kms~^ (51 pulsars). The mean veloc- 
ity of the complementary sample with ages > 10^ years is 
193 ±50 kms~^ (35 pulsars). The similarity between this 
and the overall sample is due to the great weight assigned 
to old, faint and slow moving pulsars (see also the discus- 
sion of the Iben & Tutukov results in section |5|). In Figure ^ 
we see that the distribution for the old pulsars does have a 
larger low velocity tail, as one might expect. If we place the 
velocity cutoff at 4 x 10^ years, we find 226 ± 71 kms"'^ 
for the 36 young pulsars and 198 ± 53 kms~^ for the older 
pulsars. 

Our analysis of the proper motion errors neglects the 
effects of the differential rotation or corrections to the pul- 
sar local standard of rest. However, these are not important 
corrections for two reasons. The first is that the pulsar sam- 
ple is largely restricted to distances < 3kpc, so that these 
corrections are, at most, ~ 2B x 3kpc ~ 75km s~^, where B 
is Oort's B constant. At this level it might still have observ- 
able consequences, but the second reason is that the Knax 
values are determined largely by the flux limits, rather than 
the proper motion limits. Thus, these effects are only likely 
to become important when observers attempt to extend the 
proper motion sample to greater distances. 

2.4 The Asymmetric Drift 

A reduction in the mean velocity for old pulsars has been 
noted before by several authors (Lyne & Lorimer 1994; Nice 
and Taylor 1995; Camilo, Nice & Taylor 1993). We should 
note, however, that when we refer to "old" pulsars, we refer 
to those with characteristic ages less than 10^ years, i.e. we 
don't consider millisecond pulsars because of the increased 
complexity of treating their selection effects. Nevertheless, 
in a completely model-independent way, we can demonstrate 
that these pulsars are old in a dynamical sense, because they 
show the effects of the asymmetric drift (e.g. Mihalas and 
Binney 1981). Nice and Taylor (1995) have pointed out that 
the millisecond pulsar population might possess this prop- 
erty, but it appears to be true for all pulsars with spin down 
ages > 10^ years. This is shown in Figure |^ The effect has its 
origin in the fact that any virialised population with a signifi- 
cant radial velocity dispersion will rotate about the galactic 
centre more slowly than the local circular speed (Mihalas 
and Binney 1981). 

To calculate this we restrict ourselves only to those 
pulsars with well-determined proper motions (since large 
enough error bars can reverse the sign of the transverse ve- 
locity). However, because we will not treat selection effects 
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in this case, we shall use all the proper motion pulsars that 
satisfy this and subsequent criteria, including millisecond 
pulsars. We consider a cartesian coordinate system with ori- 
gin at the sun, with positive x pointing radially outward and 
positive y pointing in the direction of £ — 270°. We consider 
the y-components of the transverse velocity, which measures 
the approximate azimuthal component of the pulsar veloc- 
ity with respect to the sun (strictly speaking, this should be 
done for each pulsar in its individual local standard of rest, 
but this rough approach demonstrates our result sufficiently 
well and remains the same if we reduce the sample radius). 
We also exclude pulsars with d > 6 kpc and within 20° of 
^=90 or 270 (where transverse velocities are primarily radial 
in the galactic frame) and of b = ±90°. This leaves us with 
37 pulsars with characteristic ages from 10^ - 10^° years. 
The signature of the asymmetric drift is thus an excess of 
positive Vy . Indeed, this is seen to striking effect in figure ^ 
where less than 10% of the pulsars with characteristic ages 
older than 4 x 10^ years have negative Vy. This indicates 
that these pulsars must be at least 10^ years old, which sup- 
ports the use of the age cutoff in section 2.3. It also provides 



a lower limit for the magnetic field/torque decay time, con- 
sistent with that found by recent analyses (Bhattacharya et 
al 1992, Hartman et al 1996). 



3 THE KICK DISTRIBUTION 

In section |^, we derived the corrected proper motion distri- 
bution appropriate to a volume limited sample. Since pulsars 
receive their kicks in a rest frame rotating about the galactic 
centre and we are interested in the low velocity tail, we need 
to consider the effect of differential galactic rotation. 

Using different initial velocity distributions, we project 
the velocities along random lines of sight in a sphere of 3 kpc 
radius about the observer and compare the inferred proper 
motion distribution with what we observe. This is shown 
in Figure ^ We find the distribution is consistent with a 
maxwellian distribution of kick velocities with velocity dis- 
persion of ~ 190 km s~^ (corresponding to a 3-D mean of 
~ 300kms~^). However, it is not consistent with the Lyne 
and Lorimer proper motion distribution or the earlier form 
suggested by Paczynski (1990)0. Neither of these analyses 
included an extensive discussion of selection effects. Thus, 
our best estimate for the kick velocity distribution is 



2 v' 



(3) 



with (Jv — 190kms~ . However, there is a significant degree 
of uncertainty about the form of the distribution, in partic- 
ular at the high end. As a demonstration of this, we propose 
a second extreme form which consists of two delta functions, 
at 250 km.s"^ and 1000 km.s"\ weighted such that 20 % of 
the pulsars acquire the higher velocity. Thus, 



p{Vk) = 0.8 5(l4-250km.s"') + 
0.2 5(14-10^ km.s"'). 



(4) 



t p(u)=4/7r(l + u^)^, where u=V/V» and V* is some normalisa- 
tion constant. We have chosen V* = 500km. s~^ to be consistent 
with the observed mean of the distribution. 



where 5{V — Vo) is the Dirac delta function. Both forms are 
shown in Figure kI 



4 LONG TERM EVOLUTION 

Our analysis above is concerned solely with those pulsars 
with characteristic ages < 10^ years. To do the same for 
the older pulsar population will require the incorporation of 
the selection effects of more pulsar surveys, as well as the 
effect of the death line and recycling of pulsars in binaries. 
Nevertheless, if we assume that all pulsars originate from 
a population with the velocity distribution (^, we may ex- 
amine the long-term evolution of this population. We have 
performed Monte Carlo simulations of such a population us- 
ing a galactic potential from Paczynski (1990). The potential 
contains two terms of the form 

(R' ^) = ,.1/2 (5) 



+ (a. + (^2H- 62)1/2)2] 



where i—1,2 represent the disk and bulge respectively. A 
third component, the halo, is represented by 



$3(r) 



GMc 
2rc 



In 1 + ■ 



+ 



2arctan (^^) 
r/rc 



(6) 



and r = (ii^ + z^^^/^ The various parameter values are 
given in Table ^. The birth positions of the pulsars are dis- 
tributed exponentially in both galactocentric radius (R) and 
disk height (2), with scale lengths of 4.5 and 0.075 kpc re- 
spectively. The integration is performed using the Burlisch- 
Stoer integration algorithm from Press et al. (1992). 

To recreate the observations we 'observe' the pulsars 
from a galactocentric radius of 8 kpc, using a volume of ra- 
dius 3 kpc about the observer (this distance is chosen to ap- 
proximately recreate the volume sampled by Harrison, Lyne 
& Anderson 1993). To make more efficient use of our simu- 
lations we can use 8 simultaneous observers spread equidis- 
tantly around the circle R = 8 kpc (see Figure bI). As long as 
the observing volumes don't overlap, the observations from 
all volumes can be added together. We calculate the trans- 
verse velocities after subtracting the circular velocity of the 
observer and derive the observed velocity distribution as a 
function of time. 

Figure ^ shows the evolution of the transverse velocity 
distribution with time. We can compare it to our corrected 
observed distributions, both for t < 10^ years. We see that 
the 10^ year distribution confirms our analysis of section ^ 
and similarly for t > 10* years (see below). 

An additional constraint on the evolved distribution 
is the magnitude of the asymmetric drift. For the kick 
distribution given by equation ^ the magnitude is small 
(Ve < lOkms"^) to begin with, reaching a maximum of 
~ 80 kms-i after 10® years and then reaching an asymp- 
totic value of ~ 60 kms"'^ after 10® years. The evolution of 
this quantity is shown in Figure ^ Again, we see we have 
the correct order of magnitude. Figure ^ demonstrates the 
magnitude of the effect (although it will be modified by the 
approach to the death line and binarity of the older pulsars) 
in our simulations as we show Vx (radial in observer frame) 
and (transverse) at different ages. Apart from the con- 
spicuous appearance of the asymmetric drift at later times. 
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Table 1. Parameters of our Galactic Model 

a b Tc M 

(kpc) (kpc) (kpc) (M0) 

Bulge 0.277 1.12 xlO^^ 

Disk 3.7 0.20 8.07 xlO^^ 

Halo 6.0 5.0 



we note also an excess of pulsars moving radially outwards 
at ages~ 10^ years. This is because more pulsars are born 
at smaller galactocentric radii and the timescale for the es- 
caping pulsars to travel a few kpc at ~ 300 kms~^ is a few 
X 10^ years. The smaller signal makes this a more difficult 
effect to observe in the real pulsar population. 

As we have noted before, a proper comparison of our 
results with the older pulsars requires an analysis of the mil- 
lisecond pulsar distribution. A treatment of the evolution of 
the entire pulsar sample is beyond the scope of this paper, 
but we can answer the question whether the oldest pulsars 
have a velocity distribution that is a direct result of the 
kick distribution given by equation To do this we model 
the effects of the surveys undertaken at Molonglo (D'Amico 
et al. 1988), JodreU Bank (Biggs & Lyne 1992), Arecibo 
(Wolszczan 1991; Nice, Taylor & Fruchter 1993; Camilo et 
al. 1993; Foster, Wolszczan & Camilo 1993; Thorsett et al. 
1993) and Parkes (Johnston et al. 1993; Bailes et al. 1994; 
Lorimer et al. 1995), although the Knax was dominated by 
the surveys from the latter two telescopes. By restricting 
ourselves to pulsars with characteristic ages > 10^ years, 
we hope to determine the velocity distribution of the old- 
est pulsars. Note that ~ 40% of this sample are still long 
period (P > 10 ms) pulsars. However, if we further restrict 
our sample to only millisecond pulsars, we obtain a similar 
result. Figure ^ shows the Vmax weighted velocity distribu- 
tion for this sample. We note that our derived distribution 
is strongly peaked between 60-100 kms~^, with the biggest 
weights attributed to the pulsars B0655-I-64, B 1952-1-29 and 
J2322-I-2057, although 6 other pulsars also lie within this 
range as well. 

A striking feature of Figure |l^ is that the old pulsars do 
not have a velocity distribution that comes from an unadul- 
terated evolution of the initial velocity distribution. How- 
ever this is not surprising, since many of the older pulsars 
are in binaries. Thus, the distribution of centre of mass ve- 
locities will be different from the kick velocities given to 
the neutron star component. Calculations of this effect have 
been performed by numerous authors for difi'erent kick dis- 
tributions (e.g. Dewey and Cordes 1986; Brandt and Pod- 
siadlowski 1995; Kalogera 1996) and introduce further un- 
certainties into the problem through the initial distribution 
in orbital period, companion masses and pre-supernova he- 
lium star masses. Thus we shall defer detailed discussions 
to future work, although we note that the effect is indeed 
to reduce the high velocity tail. Another effect is that the 
low velocity end is also expected to be depleted due to the 
gravitational scattering off giant molecular clouds (Spitzer & 
Schwarzschild 1951, 1953) or spiral arms (Barbanis & Wolt- 
jer 1967). This is a concern for calculations of the observ- 
ability of old neutron stars accreting from the ISM ( because 
the accretion rate and thus luminosity oc v~^), and Madau 
& Blaes (1994) have recently calculated this efi'ect, yielding a 



peak in the range 60-90 kms ^ after 10^ years ( Although 
they used the Narayan & Ostriker (1990) kick distribution). 



5 DISCUSSION 

Our approach above is designed to provide a robust esti- 
mate of the characteristic birth velocities of pulsars. We 
have tested the procedure with simple Monte Carlo simu- 
lations using known distribution functions. We find that we 
can reproduce the mean velocities to within ~ 30 — 50 km 
s~^, although the exact shape of the distribution at the high 
velocity end (> 300kms~^) is not well constrained. To bet- 
ter estimate the shape of the distribution will require more 
detailed modelling and the use of more of the pulsar popula- 
tion information, such as was done by Narayan and Ostriker 
(1991) or Bhattacharya et al. (1991). 

Another recent analysis by Iben and Tutukov (1996) 
finds a large birthrate of very slow (< lOkms^^) pulsars. 
Their treatment neglects flux limits (although some infer- 
ences are made on the basis of a nearby sample only, their 
full analysis uses pulsars at all distances), and treated the 
proper motion limits using various ad hoc analytic cutoffs. 
They also used no upper age cutoffs, so that their large 
birthrate of slow pulsars was due to 3 old (t > 10^ years) 
pulsars, which acquired significant weight because of their 
small luminosities. When restricting ourselves to ages < 10^ 
years, we do not find any such low velocity tail. Furthermore, 
we find that the flux cutoff is responsible for the detection 
limit of more pulsars than the proper motion limit. 

Some constraints on the nature of pulsar kicks can be 
obtained by analysing the properties of binaries containing 
neutron stars. Recently, in the light of the results of Lorimer 
and Lyne (1994), Brandt and Podsiadlowski (1995) anal- 
ysed the effect of the revised kick distribution on the post- 
supernova orbital parameters of neutron star binaries. Their 
analysis indicated that an isotropically distributed kick ve- 
locity of 450 km s~^ was inconsistent with the eccentricity- 
orbital period distribution of the observed binary popula- 
tion. However, a kick velocity of 200 km s~^ was perfectly 
consistent. Similarly, Wijers et al. (1992) obtained an upper 
limit of 400 km s~^ for a characteristic kick velocity from 
an analysis of the eccentricities of the known double neu- 
tron star binaries. While these results are in good agreement 
with our analysis, we should note that Brandt & Podsiad- 
lowski showed that the Lyne & Lorimer distribution, despite 
it's high characteristic velocity, was also consistent with the 
period-eccentricity relation, by virtue of it's large low ve- 
locity tail. Thus, the surviving binary parameters do not 
provide a means for distinguishing between the two distri- 
butions. 

The most direct test of the pulsar kick velocities is to 
flnd the supernova remnant associated with the birth of a 
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given pulsar. If the pulsar was born at the centre of a circu- 
larly symmetric supernova remnant, we may infer a proper 
motion from it's displacement and thus obtain a velocity. 
Much work has been devoted to this (Caraveo 1993; Frail, 
Goss & Whiteoak 1994), and the results indicate velocities 
significantly larger than the mean velocity derived by other 
methods. In fact, the mean velocity of the Frail et al. sam- 
ple is 990 km s~^, with a median of 480 km s~^. Whether 
or not this discrepancy with respect to our results is real 
depends on the veracity of the various assumptions used 
to infer a proper motion from a supernova association (see 
Kaspi 1996). Some of the problematic assumptions discussed 
by Frail et al. include the difficulty of defining the shape 
of a remnant and the possible displacement of the super- 
nova blast centre with respect to the geometric centre of 
the remnant (perhaps caused by expansion into an inhomo- 
geneous surrounding medium). It is illuminating, although 
not conclusive, to note that the two pulsars associated with 
supernova remnants which have measured proper motions 
are the Crab and Vela pulsars, with transverse velocities of 
150 and 120 km s~^ respectively (In particular. Frail et al. 
note that the methods used on other pulsar-remnant asso- 
ciations would imply a velocity of 800 km s^^ for the Vela 
pulsar!) 

A revised velocity distribution can possibly affect the 
results of statistical analyses of the pulsar population as 
a whole. Indeed, it may even contribute to an explanation 
for the conflicting claims concerning magnetic field decay. 
Narayan and Ostriker (1990) used a maxwellian distribu- 
tion for their paper on the pulsar population as a whole. 
They used two different populations of pulsars and their 
kick distribution was a function of magnetic field. Their two 
populations had velocities that varied from 80-250 kms^^ 
for their S population and from 20-100 kms~^ for their F 
population. The analysis of Bhattacharya et al. (1991) also 
used a maxwellian, but with a somewhat lower dispersion of 
110 km s~^. Our value lies a little above these characteristic 
values, consistent with the fact that the dispersion measure 
distances are now thought to be larger than those used in the 
above analyses. More recent analyses, such as that of Hart- 
man (1996), attempt to constrain the velocity distribution 
itself using this approach. Results from such studies may dif- 
fer significantly from ours because of the sundry additional 
assumptions required. In particular, the Crab and Vela pul- 
sars acquire significant weight in pulsar current studies (see 
Phinney & Blandford 1981) because they require significant 
luminosity evolution to avoid a large population of high field 
objects near the death line. This could produce a larger low 
velocity tail than in our case. Similarly, inclusion of some of 
the more optomistic pulsar-supernova remnant associations 
in the population synthesis can result in a high velocity tail. 
Relative to such analyses, our results represent a minimal 
model, required to describe the observed proper motions un- 
der the assumption that the death line does not remove a 
significant fraction of pulsars before a spin-down age of 10^ 
years. 

An interesting feature of our result is the lack of a pro- 
nounced low velocity tail. This has implications for the re- 
tention fraction of neutron stars in globular clusters, as well 
as the existence of some obviously low velocity pulsars with 
ages > lO'^ years. Globular cluster central escape velocities 
are < 50 km s~^. If all kick velocities are ~ 250 km s~^, then 



no pulsars born from isolated stars are retained! If the dis- 
tribution is maxwellian, then the fraction retained is about 
0.2%, which is still extremely low! Yet, the retention frac- 
tion of neutron stars is claimed to be of the order of 10% 
(Phinney 1993) or higher (Hut and Verbunt 1983). However, 
it is possible that binaries, either primordial or dynamically 
formed, could be responsible for the pulsars found in globu- 
lar clusters (see Hut et al. 1992; Drukier 1996). Brandt and 
Podsiadlowski (1995) find that ~ 17% of binaries remain 
bound if the kick velocity is 200 km s~^, which is the right 
order of magnitude to explain the required mass in dark mas- 
sive remnants. However, this is likely to be an upper limit 
on the retention fraction because even systems which re- 
main bound can receive significant centre-of-mass velocities. 
Thus, Drukier (1996), using the results of Lyne & Lorimer 
and Brandt & Podsiadlowski, has shown that most globular 
clusters would retain ~ 1 — 5% of their neutron stars. In the 
light of this, we should point out that the lack of a low ve- 
locity tail in our distribution is not an artefact of the Knax 
weighting. Of the 51 pulsars in our sample with ages less 
than 10^ years, the lowest transverse velocity is 70 km s~^. 
If there are young pulsars with very small velocities, then 
they have not been measured yet. Another possible compli- 
cation is the creation of fast (P < 0.1 s) pulsars with initial 
timing ages > lO'^ years. The birthrate of such pulsars is not 
constrained by our analysis because of both our age cutoff 
and the restriction of our analysis to the early Molonglo and 
Green Bank surveys. 

In conclusion, we have shown that the distribution of 
young pulsar proper motions, corrected for selection effects, 
is consistent with a characteristic kick velocity at birth of ~ 
250-300 km s"^ We find little evidence for a significant low 
velocity tail to the kick distribution. Our method is robust 
in its reproduction of the mean velocity and low velocity 
shape of the distribution. However, the shape of the distri- 
bution at velocities > 300 km s~^ is not well constrained 
by this method because of poor statistics. Our results are 
in good agreement with the properties of binaries contain- 
ing neutron stars and pleasantly close to the value inferred 
from numerical supernova simulations (e.g. Burrows, Hayes 
& Fryxell 1996, who inferred kicks of ~ 300kms~^ from 
core recoil during the collapse, but neglected asymmetries 
in the initial mass distribution, see Burrows & Hayes 1996, 
which could lead to even higher velocities) . Finally, we have 
shown that the pulsars with long characteristic ages show 
the asymmetric drift, corresponding to a dynamically old 
population. The velocity distribution of these pulsars is af- 
fected by their genesis in binaries as well as their subsequent 
motion through the galaxy. 
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Figure 1. Sample Definition: The filled circles indicate the 
pulsars that we include in our analysis. The open circles de- 
note those that are excluded. Since the purpose here is simply 
to demonstrate which pulsars are in the sample discussed, we 
omit any error bars. Of the three excluded pulsars with P > 0.1s, 
two are recycled binaries, 30655+64 and B0820+02. 



Figure 2. Luminosity and Velocity for the Proper Mo- 
tion Pulsars: We again include all the pulsars with proper mo- 
tions on this plot. The filled circles will be the ones to which our 
analysis applies. The dotted line indicates a proper motion of 5 
mas yr~^ and fiux 4 mjy. This line does not represent a cutoff 
over most of this diagram because making a pulsar brighter at 
a given distance will move it to the right and making a pulsar 
faster at a given distance will move it up, thus one can popu- 
late both sides of the line with observable pulsars. Nevertheless, 
at the high luminosity/high velocity end, it should represent the 
limiting case. It is also noticeable that, on average, higher velocity 
pulsars have higher luminosities, and so will be overrepresented 
in an unweighted sample. 



Figure 3. Selection Effects: The solid circles represent those 
pulsars detected in the Green Bank surveys. The open circles are 
the ones detected in the Molonglo surveys. The dotted lines indi- 
cate the approximate limiting fiux as a function of P for DAf =50 
cm~^ pc and the most sensitive Molonglo and Green Bank sur- 
veys. The dashed lines are for DM=200 cm~^ pc. 



Figure 4. Detectable Volume for Each Pulsar: Pulsars with 
small luminosities will only be observable near the galactic plane, 
so that their Vmax will be spherical (neglecting other selection 
effects for the moment). Pulsars with large luminosities will be 
observable much further away, out to the limits of the disk that 
such a population, born in the galactic plane, would fill. In this 
case, the spherical Vmax will be cut off above the limits of the 
disk height given hy z = Vz x t. 



Figure 5. The Weighted Proper Motion Distribution: The 

size of each circle is proportional to the logarithm of the weight 
accorded that pulsar. The vertical dotted line indicates the di- 
viding line between what we consider 'young' and 'old' pulsars. 
The horizontal dashed lines indicate the range of values we in- 
fer for the mean proper motion of the pulsars in each sample. 
The short horizontal solid line at left indicates the mean proper 
motion value quoted by Lyne & Lorimer (1994). 



Figure 6. The Inferred Proper Motion Distribution: The 

cumulative probability distribution P(y < Vt) is shown for difi'er- 
ent degrees of correction. The dotted line shows the entire sample 
with no Vmax weighting and no treatment of upper bounds. The 
short dashed line is for only the 66 pulsars with detected proper 
motions corrected using our Vmax correction (i.e., all those with 
only upper limits were left out) . The long dashed line is the prop- 
erly corrected sample (with V^ax corrections and Kaplan-Maier 
estimator) but only those with characteristic ages greater than 
lO'' years. The solid line is the complementary sample of only 
those pulsars with characteristic ages less than lO'^ years. The 
95% confidence levels for the corrected distributions lead to an 
uncertainty of ~ 0.05 in P for velocities < 300 km s~^. Above 
that, the statistics become uncertain and we cannot say much 
about the distribution. 



Figure 7. The Asymmetric Drift: The solid squares indicate 
a positive Vy and open squares indicate a negative Vy . The circled 
points indicate binaries. The horizontal dashed line indicates the 
order of magnitude of the sun's motion within its local standard 
of rest, and velocities below this will experience a contamination 
of the asymmetric drift. The vertical dashed line is at 4 xlO® 
years, and represents an approximate division between "young" 
and "old" pulsars, namely those that show the asymmetric drift 
and those that don't. Of the young pulsars, only 33% have positive 
Vy, while 91% of the old pulsars have positive Vy. The dashed 
line indicates the evolution of the mean asymmetric drift velocity 
obtained from the calculation in section hi 



Figure 8. Determining the Birth Velocity Distribution: 

The heavy solid line is the corrected proper motion distribution, 
with the formal uncertainty shown by the shaded region. The 
short dashed line is the proper motion distribution obtained for 
a distribution given by two delta functions at 250 km/s and 1000 
km/s, in random directions. This distribution and the maxwellian 
given by (^ (thin solid line) are consistent with the data. The 
alternating and long dashed lines represent the results of Paczyn- 
ski (1990) and Lyne & Lorimer (1994), neither of which are fully 
consistent with the data. 

Figure 9. Monte Carlo Pulsars: The right half of this pro- 
jection onto the galactic plane is for the initial positions of the 
pulsars. The mirror images of the positions of the same pulsars 
after 10^ years is shown on the left hand side. The circles indicate 
the volumes sampled by the 'observations'. 

Figure 10. Velocity Distribution Evolution: The solid lines, 
labelled by log(age), are the time-dependant transverse velocity 
distributions produced by the initial three dimensional distribu- 
tion labelled Vq (dotted line). The dashed line indicates the ob- 
served 'young' (< lO'^ years) pulsar proper motion distribution. 
The heavy line indicates the observed 'old' (> 10* years) proper 
motion distribution. The line is solid where the contribution is 
from a binary and dotted when the corresponding pulsar is iso- 
lated. 



Figure 11. The four panels (all to the same scale) show the 
radial (Vx) and transverse (Vy) velocities in the observer frame 
at four different ages. After 10^ years the velocity distribution 
is isotropic because only locally born pulsars are seen. After lO'^ 
years an excess of radially outward motion is seen due to the net 
fiux of pulsars moving out from the inner parts of the galaxy. 
After 10* years only pulsars from the lower velocity tail are seen 
and these clearly exhibit the asymmetric drift. 



© 1997 RAS, MNRAS 000, 0-?? 





1 


1 1 1 1 1 M 


1 III 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1114 

• 


1000 












o 






• ••it* * 


1 

m 






o 




B 100 




°oo 

o 


• 


> 






o 




• 

. • • • 

• • 


10 


= o 

1 





1 III 


• 

1 1 1 1 



0.001 0.01 0.1 

P (s) 



- 5 - 




rm \ \ I I I 1 1 1 1 \ \ I I I 1 1 1 1 \ \ I I I 1 1 1 1 \ \ I I I 1 1 1 1 \ \ r 



1000 



CO 

a 

> 



® 



100 - 



10 



© 



C) 



C) 



C) 

® c: 



C) 



O 



O 



o 



-U I I I I I I I 



J I I I I I 1 1 



o 



o 



T 



O 



o 




CD 



O 



() 



I M I I 



c 



e 



b 

C) 



o 



(3) 



J l_L 



C) 



o 



o 



o 



11 



o 







o 



J I 1= 



104 



10^ 10^ 
t (years) 



10^ 



1 ^ ^ l^r'-i^^ 



1 777T. 



0.8 



0.6 



> 

V 

-I- 



0.4 



J 

I r 



0.2 







— r 

,f---' : 

_^JJ ^ : ■ 

_L1 I I I I I l_L 



No correction 

V„>ax. No KM . 
T>10^ years 
T<10^ years 



J I I I I \ L_L 



J L 



100 



V (km.s-1) 



1000 



"1 I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I r 



1000 



: 5/10 



[] 



i 



[] 

[] 



100 



"[] 



20/2 



i4 



1957+20 



1855+09 



[] 



10 



1937+21 

J I I I I I I I I I I I I I I I I I I I I I I I I I I I L 



105 



106 



10^ 108 

T(years) 



109 



10 



10 



20 



"1 I I r 



"1 I r 
o 



"1 I I r 



"1 I r 



r 



t = 



o 



10 







-10 



eg 



o o 
o 



o 
o 



o .-■ 



I • - • ' • - 



o " o a o I* •W 



8 



CP o,o-, 



O' 



o 
o 



o o . 

O ', fP 



oo 
o 'o 



q.- © 
o ^ 



>•••: /•* * • • 



-.On 

° O \ O 
O O 



• • ••• I 

•t • • 



O, 



O • 
O'. 
O ' 
O 



-S1 



20 



t = 109 



J I I L 



O 
o 

J I L 



o 
o 



J I I L 



J I I L 



-20 



-10 





X (kpc) 



10 



20 







100 

(km.s-1) 



1000 



